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I. INTRODUCTION 

The study of systems composed of a large number of 
degrees of freedom was endowed a definitive means of rea- 
soning with the introduction of the differential stochastic 
dynamical framework about a century ago Specifi- 
cally, differential stochastic dynamics allows one to de- 
scribe the time evolution of a given observable as a func- 
tion of macroscopic variables of the system as well as 
non-deterministic ones. The latter are expressed in terms 
of noise(s), reflecting the microscopic features of the sys- 
tem, whose univocal description is beyond the bounds of 
possibility 0, Q. Despite the impossibility of a deter- 
ministic description of the evolution of the observable, 
one can successfully obtain time dependent statistical 
details. Undeniably, Einstein's and Bachelier's groun- 
breaking works, respectively on Brownian motion [J] and 
stock price movements are outstanding examples of 
the relevance of differential stochastic dynamics in di- 
verse fields. 

The inference of a differential stochastic process is 
generically made from statistical features such as statis- 
tical moments, correlation functions and probability den- 
sity functions built from a time series of measurements. 
In particular, a well established theoretical background 
for identifying the stochastic dynamics exists for the im- 
portant class of processes following the Markov prop- 
erty [HQ- However, in practice, several hindrances arise. 
Specifically, besides matters related to the finite size of 
data sets [6| and to the direct error of the measurements 
(associated with the quality of the equipments and/or 
their calibration Q ) , the sample rating (the spell between 
logged measurements) plays a crucial role in the deter- 
mination of actual underlying stochastic process. As a 
matter of fact, we are generally left a set of snapshots 
reproducing a fraction of the events that occurred in the 
continuous time process from which the dynamics is due 
to be determined. In any case, it would be important to 



compare the acquisition interval with the characteristic 
times of the process to determine whether the Kramers- 
Moyal (KM) coefficients estimates are trustworthy. In 
this respect, another difficulty is that the characteristic 
times are not always known beforehand, although they 
can be estimated, e.g., through the computation of auto- 
correlation functions. Moreover, when the sampling in- 
terval is found to be too long compared with the char- 
acteristic timescales, for suitably uncovering the process, 
it is not always possible to upgrade it, particularly for 
historical data. For such cases, it is therefore of primary 
importance a careful analysis of the impact of the time 
interval r of data sampling on the observable finite-time 
KM coefficients, specially as one approaches the indepen- 
dence limit. 

For a Markovian timeseries, one can obtain the evo- 
lution equation for the conditioned probability density 
(PDF) by computing the KM coefficients, 

D k (x ) = lim D k (x ,T) , (1) 

with 

D k {x ,T) = -j- J <IxP(x,t\xo,0)[x(t) -x ] k 

= ^t{ k j )(^r\ (2) 

where we have already assumed that the coefficients are 
(at least locally) stationary (for non-stationary data-sets 
see Ref. @). Aiming to simplify the notation, we denote 
the statistical averages conditioned to the initial value 
xq: (• • • )|a;=!Eo = ('■'); w hd e we will reserve (•••)« for 
usual stationary unconditioned averages. 

In practice, only the finite-time estimates Dk(xo,r) 
can be directly computed, with r limited by the minimal 
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time interval, T m i n , of data acquisition. Furthermore, 
when the sampling interval t is larger than the correla- 
tion time, i.e., the PDF P(x,t\xo, 0) becomes uncondi- 
tioned and under stationarity, one gets, 



j^indep 



(xq,t) 



^- J dxP(x)[x - a*]' 



fc!< 



: ([x(t) -X ] k )u, 



(3) 



Notice that, in Eq. @ the stationary averages are un- 
conditioned, hence do not depend on x$, differently from 
Eq. ([2J). Therefore, Eq. (|3]) represents a fc-ordcr poly- 
nomial in xq, no matter how complex the intrinsic co- 
efficients Dk are. In practical applications, this feature 
introduces a complete uncertainty on the actual form of 
the intrinsic coefficients, in the absence of a model a pri- 
ori m. 

In this manuscript we obtain the finite-time estimates 
for several standard processes described by Ito-Langevin 
stochastic differential equations, 



II. ESTIMATING CONDITIONAL MOMENTS 

Given the Fokker-Planck equation (FPE) for P = 
P(x,t\x o ,0) 



d t P = -d x [D 1 P] + d xx [D 2 P] 



(5) 



the evolution equation for a mean value can be obtained 
by multiplying both sides of Eq. (J5J) times the quantity 
to be averaged and integrating by parts with suitable 
boundary conditions (vanishing at the boundaries). In 
this way, the equations for the moments have the form 



d(x n ) 
dt 



n(x n - 1 D 1 (x)) + n(n - l)(x n ~ 2 D 2 {x)) . (6) 



If Di(x) is linear and D 2 (x) is at most quadratic, the 
equations for the moments can be successively solved 
from the lowest order n = 1 [Toj . However, a hierarchy of 
equations depending on higher order moments generally 
arises. In such cases, one can still resort to approximate 
techniques such as hierarchy truncation, substitution of 
terms, etc. Alternatively, the conditional averages can 
be computed by means of the short-time solution expan- 
sion 0, which leads to 

(Q(x))(x ,t) = ^2\tf(x)] n Q(x)\ Xo r n /n\ , (7) 



dx = D x {x)dt + y/2D 2 (x)dW t 



(4) 



where Wt is a standard Wiener process. In previous 
work Q, this task was accomplished for the particular 
class of processes with D\ = —ax and D 2 = Ax 2 + C, 
by means of Ito- Taylor expansions. Herein, we will con- 
sider a larger class of processes following two different 
approaches. The first one, based on the solution of the 
evolution equation of the raw statistical moments, rep- 
resents a simplification regarding the whole process of 
appraising the finite-time coefficients when the drift is 
linear. Second, we employ the Fokker-Planck adjoint op- 
erator technique that is useful for processes with nonlin- 
ear drifts, for instance. 

Furthermore, we will analyze the extreme instances 
of no-fluctuations and independence. The deterministic 
limit (absence of fluctuations) is relevant in stochastic 
dynamical systems since it provides a reference on the 
functional behavior of quantities that can be perturbed 
by augmenting the intensity (broadness) of the noise. In 
the meanwhile, the independence limit sets the bounds 
expected for large acquisition interval and/or large noise 
intensity. 

Although it is not always possible a direct inversion 
of the finite-time estimates to extract the intrinsic ones, 
the expressions herein obtained can be contrasted against 
the observed KM coefficients, allowing the identification 
of the underlying process as well as the associated pa- 
rameter values. 



where L'(x) = D\d x +D 2 d xx is the (backwards) Fokker- 
Planck adjoint operator (see also 

In the deterministic case, it is enough to solve Eq. © 
for the first moment, which becomes 



d(x) 



det 



dt 



(8) 



where the superindex stands for "deterministic" and 
whence we obtain the full range of values, 



lt k\det 



{(x) ) . Then, the deterministic part of Dk, with k > 1, 
is 



f)det 



I 



det 



r fe-l 



-(Df et ) k 



(9) 



Let us note that the deterministic Dk are non-null func- 
tions, although they vanish in the limit r — > except for 
k = 1. Nevertheless, we should bear in mind that for 
not too small r the deterministic part can dominate over 
the noise contribution. This represents another draw- 
back one may face in a practical application. In fact, if 
the noise is small, the deterministic contribution to Dk 
may screen the information on the noise. 



III. LINEAR DRIFT 

Let us first consider the case of linear drift, 

D 1 (x) = -ax + b, (10) 
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with a > assuring the existence of a stationary solution. 
As a matter of fact, a raft of phenomena are described 
by an exponential relaxation, which is concomitant with 
a parabolic drift potential [|| . From Eq. © one has 



d(x) 



= —a(x) + b , 



whose solution is 



(x(t)) = x z + -(1 - z) , 
a 



(11) 



(12) 



assuming the initial condition (x(0)) = xq and defining 
z = c~ aT . Hence, according to Eq. ©, 



Di(x ,t) = - I x ) - — 

a J t 



(13) 



that in the limit t — > recovers D\. The opposite limit 
of independence yields — (xq — b/aj/r in accordance with 
Eq. ©. 

It is noteworthy that Eq. (IT31 is independent of the 
particular form of D 2 . This is due to the fact that in 
(and only in) the linear case, Eq. (ITTj) coincides with its 
deterministic version, Eq. flSJ), hence D 1 = Df et . For 
other forms of the deterministic part of the stochastic 
differential equation, the derivative of (x) with respect 
to the time is equal to more complex expressions which 
depend on higher-order moments and/or powers of the 
first conditional moment. Despite the slope evolution 
with r, it is worth noting that the finite-r drift preserves 
the linear form of the intrinsic D\. 

Equation (|12[) also allows one to estimate the condi- 
tional two-time covariance function given by Q, 

K(t,xq) = J dx xxo P(x, t\xo, 0) = xq(x(t)) 

= xlz + x -(l-z), (14) 
a 

that exponentially decays towards the long-time average. 
Still one has to average over initial conditions xq to find, 

K{t) = J dxdxo xxqP(x, t; Xo, 0) = J dxoK(r, xq)P(xo), 

where P{xq) is usually the steady PDF. However, any 
choice of P{xq) will not alter the exponential character 
of the correlations since the integration is not carried out 
over time. In addition, the correlations are only ruled by 
the drift and do not depend on the noise amplitude. It 
is straightforward to understand that this quirk comes 
to pass because of the functional properties of the condi- 
tional average. 



1. Second and higher order conditional moments for 
quadratic noise intensity 

If together with drift linearity the noise intensity is 
second-order polynomial, 



D 2 (x) = Ax 2 +Bx + C. 



(15) 



with parameters such that D%(x) > 0, then each evolu- 
tion Eq. © is linearly coupled to those of lower order 
only. In this case the equations can be solved one after 
another. In particular, for the second moment, one has, 



d(x 2 
dt 



= 2(-{a-A)(x 2 ) + {b + B)(x)+C). (16) 



Thus, we just need to solve a non-homogeneous 1st order 
differential equation of the form 

/' = Af + Bc~ at + C, (17) 

whose solution, given the initial condition /(0), is 

= aC + A(B + C+(a + A)f(0)) &At 



A(a + A) 



B - at C 
-e 



a + A A 
Therefore, (x 2 (t)) = /(r) with the identifications 

A = -2(a-A), 
B = 2(6 + B){x -b/a), 
C = 2(b + B)b/a + 2C, 



(18) 



and a — a, following Eqs. (fT6| and (|T2|) . and the initial 
condition /(0) = x\. Namely, 



(X (T)) — * 

b(b + B) 
+ a(a - A) (a - 2 A) 



1 1 r, ^ + P I 2 \ 

WX n + 2 T-r(Z — Z W)Xq 



a-2A 

(az 2 w - 2(a - A)z + a - 2A) 



C 



A 



(1 - z 2 w) 



(19) 



— r,2AT 



where z = e~ aT and w = 

The finite-r second KM coefficient can be obtained 
by means of Eq. © with (x(t)) and (x 2 (t)) given by 
Eqs. (|T2"f and ([T^)) . respectively. Explicitly, we find 

D 2 (x , t) = — ( [1 - 2z + z 2 w] xl 



+2 



—(z-z 2 w) + -(z-l) 

a — 2A a 



xo 



+ 



b{b + B)(az 2 w - 2 (a -A)z + a- 2 A) 



i(a - A) (a - 2A) 



C 



-A 



(1 - z 2 w) 



(20) 



When t — y 0, D2 — ^ D 2 , while for ar >> 1, Eq. ([3]| is 
verified. 

These results permit one to embrace many funda- 
mental models such as Ornstcin-Uhlenbeck, Feller and 
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harmonic drift with additivc-(lincar)multiplicative noise. 
Also in finance, in most well-known models of volatility 
(a = y/x), the drift is linear and D 2 = D\x\^ 7 with j3 = 
(Ornstein-Uhlenbeck process), 1 (square-root, Feller 
or Cox-Ingersoll-Ross model) and 2 (Hull and White 
model). 

Let us also recall that for the linear drift, its finite- 
time expressions are also linear in xq. Similarly, with 
linear drift and /3 = 0, 1, 2 the dependence of D 2 on xo is 
quadratic. As a matter of fact, a similar scenario holds 
for fci/i-ordcr finite-time KM coefficients. In other words, 
the kth-oidei KM coefficients are polynomials of order 
k, whose coefficients depend on r. Thus, it is not simple 
to nimbly rescue the value of /3 just by identifying the 
polynomial order of the observed coefficients. 

For the quadratic noise intensity, the evolution equa- 
tions of higher order conditional moments, from Eq. ([6]), 
are 



d(x n ) 



dt 



A n (x n ) + B n (x"- 1 ) + C n (x n - 2 ) , (21) 



where 



4(34- a), B = l2C 2 /(A-a) + l2Cx 2 , C = -12C 2 /(A- 
a) and a = 2 (a — A), together with the initial condition 
/(0) = x o, yielding 

, a. .. a r a 6C(z 2 w — z 4 w 6 ) , 



5A 



3C" 



(a - A)z 4 w 6 - 2(a - 3A)z 2 w + a - 5A 
(a- A)(a - 3A)(a - 5A) 



(24) 



Thus, 



= 2^1 (^ 4 (r))-4xo(x 3 (r)) 



+6x 2 {x 2 (t))-4x 3 (x(t))+x 4 



(25) 



is obtained by substitution of the averaged quantities 
Eqs. dHJ|, (ED), QU) and (JUJ). The present results gen- 
eralize those previously obtained for the particular case 
Dt = -ax (with a > 0) and D 2 = Ax 2 + C in (<|. 

(ii) For the Feller process (A = C — 0), by successive 
integration of Eqs. (|2ip and using (J5|) we find, 



A n 



n [(n - 1) A - a] ; 
n [(n- 1)5 + 6]; 



= n (n - 1) C. 



High-order KM coefficients are particularly relevant as 
well. In other words, for the present class of Marko- 
vian processes the Pawula theorem Q implies that they 
must vanish. In the absence of any other reason to dis- 
card the validity of Chapman-Kolmogorov approach, if 
the finite-time conditional higher-moments are non-null, 
it is important to probe whether it is the outcome of a 
finite-time effect. Let us restrain our calculations to the 
explicit formulae of the third and fourth coefficients in 
two important cases: (i) the reversion to a null value 
in the absence of coupling between additive and (linear) 
multiplicative white noises (6 = and B = 0) and (ii) 
the Feller process. 

(i) For the additive-multiplicative process with 6 = 
and B = 0, notice that Eq. © with n = 3 has the 
same form of Eq. (fT7|) . through the identification A = 
— 3(a — 2 A), B = 6Cxq, C = and a = a, and the initial 
condition /(0) = Xq, yielding 



<z 3 (r)> = ^j{z ~ {zwf) + xl{zwf . (22) 
a — SA 

Then, according to Eq. ([2]), one has 
1 



D 3 (x ,t) = —^{x 6 (T)}-3x (x'(r)} + 3x' (x(T))-x 6 j , 

(23) 



where the averaged quantities are given by Eqs. ([22]). (JT9J) 
and (fI"2"T) . respectively. 

For the fourth moment, Eq. ^ with n — 4 has the 
same form of Eq. (fTTj) . through the identification A — 



D 3 (x ,t) 
-3 
and 



a-^) 3 ( x 3 1 b-(b + 2B)z 
(b + B)(b- (b + 2B)z) b(b + B)(b + 2B) 



a 2 (l-z) 



Da(x ,t) 



24r 



-X 



xt-A 



6- (b + 3B)z 



(26) 



a(l — z) 



+6 



(6 + 2B)\(b + 3B)z - 2(6 + B)]z + 6(6 + B) 
a 2 (l-z) 2 ; 



A (b + B){b + 2B)(b- (b + 3B)z) 

3(1 - Z) X ° 

b(b + B){b + 2B)(b + 3B) 



(27) 



Considering that for the Feller process (x)™ = 
(B/A) n r< ^(b/~B) > ' then the KM coefficients can be rewrit- 
ten in terms of the unconditioned raw moments, which 
can be straightforwardly estimated from data, as follows, 



~ (1 — z) , 
Di{x 0i t) = {x- x ) l 



D 2 (x ,t)= { ±^((x-x ) 2 ) u , 



D 3 (x , t) = — ^ Ux - X f 



Gr 



2 ~2 



z(l-z) a 



(x(x - X )) u Xq, 
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FIG. 2: Color online. Impact of the noise intensity on the 
first and second finite-time KM coefficients for the same pro- 
cess of Fig. [1] In all cases r = 5 and the different values of 
B are indicated on the figure. Black solid lines correspond 
to the intrinsic coefficients, black dashed lines to the inde- 
pendence limit forms at r = 5 and black dotted lines to the 
deterministic forms (Eq. ([9])) at r = 5. Symbols correspond 
to numerical estimates and the associated colored (thin) solid 
lines to the theoretical expressions for finite-time coefficients 
given by Eqs. ()13|) an d (|20p . In the lower panel the inset 
shows the same data in log-linear scales. 



FIG. 1: Color online. Impact of r on first, second and fourth 
finite-time KM coefficients for the process dx = (— ax + b)dt + 
V2Bx dW, with a = 0.1, b = 1.0 and B = 0.25. The values of 
r used are indicated on the figure. Artificial timeseries were 
generated by means of the Euler algorithm with dt — 10~ 3 
and 10 6 data points of the timeseries were considered for the 
computation of KM coefficients in each case. Symbols corre- 
spond to numerical estimates and the associated colored solid 
lines to the theoretical expressions for finite-time coefficients 
given by Eqs. ((T3j) . ((20)) and ((27)) . Black solid lines corre- 
spond to the intrinsic coefficients D\ = — ax + b, D2 = Bx 
and D4 = and black dashed lines (practically coinciding 
with the curves for r = 200), to the independence limit forms 
given by Eq. © with r = 200 



and 

Di(x ,T) = ^ 94 ^ ((x- xof) u 

z(l-z) 3 <jI , 2 \ z 2 (l-z) 2 ai 2 

2r W„ {x)u 



where a 2 x = (x 2 ) u - {xf u . 

In Fig. [TJ we exemplify the behavior of D\ , D2 and D4 
for the Feller process dx = (—ax + b) dt+v2Bx dW (that 
is A = C = 0, or also f3 = 1) as the independence limit is 
approached. Notice the variety of behaviors that can be 
observed even for not too large r. The slope of the linear 
D\ changes with r taking the values — (1— z)/t, while L>2 
soon becomes quadratic with t. Also D4 largely departs 
from the intrinsic null value. For all the coefficients, the 
independence limit is practically attained at r ~ 200 
(r » l/a = 10). 

Figure [5] illustrates the impact of the intensity of noise 
on the finite-time estimates. Noise intensity does not af- 
fect D\ which coincides with the deterministic straight 
line. D2 largely disagrees with the true form D2, de- 
parting towards the independence limit, for large noise 
intensity, and towards the deterministic limit for small 
noise. Of course, as one approaches the deterministic 
(no-fluctuations) limit the range of xq shrinks around the 
deterministic (equilibrium) value. 
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When the drift is linear, whatever the noise intensity, 
the parameters a and b can be unraveled by a linear fit to 
the observed D\ : the value of the slope — (1 — z)/t allows 
one to determine the parameter a, once r is known, while 
the abscissa at which D\ vanishes corresponds to b/a. 
That is, if the observed drift can be considered linear in 
good approximation, a direct inversion of D\ to obtain 
the intrinsic D\ is possible. This will also facilitate the 
obtcntion of the functional form (and parameter values) 
of D 2 . 



2. Non-quadratic noise intensity 

For f3 ^ 0,1,2, or more generally for non-quadratic 
D2, the finite-time D 2 are also non quadratic in x. The 
equation for the second moment, following Eq. ©, 

^p- = -2a{x 2 ) + 2b(x) + 2(D 2 (x)) (28) 

cannot be solved straightforwardly in this case. 

Let us consider the especial case with /? = 3 (known 
as 3/2-model) used as a volatility model as well [12(- For 
small noise intensity, one can find the correction c 2 to 
the deterministic part of {x 2 } by means of expansion ([7|. 
taking Q(x) = x . By identifying the general form of the 
coefficients of r" , with the aid of algebraic manipulation 
programs, we obtain 

C 

c 2 = — (b 3 + 6b 2 {ax -b)z+{3b 3 -6a 2 bxl + 2a 3 x 3 l )z 2 
+ 6abT(ax - bfz 2 - 2z 3 {ax - b) 3 ) . (29) 

Then 

(x 2 ) = ((x) det ) 2 + c 2 + 0(C 2 ) (30) 
and on that account, 

D2 = l(D d 1 et ) 2 + ^ + 0(C 2 ), (31) 

which in the limit r — > tends to D 2 = Cx 3 as expected. 



IV. NONLINEAR DRIFT 

As in the case of linear drift and non-quadratic noise, if 
the drift is nonlinear, the equations for the moments can- 
not be solved successively from the lowest order. More- 
over, the useful independence of the evolution equation 
for (x (t)) on higher order statistical moments, which are 
noise dependent, fails. Consequently, the existence of a 
nonlinear drift introduces a troublesome relation between 
E>i and the noise as well as a dependence of the corre- 
lation function on the same noise. However, we can still 
use the calculations at the deterministic limit to repre- 
sent the upper bound which is quite reliable for the cases 



presenting small noise intensity. Explicitly, if the contri- 
bution of diffusion is neglected, the evolution is almost 
deterministic. Accordingly, we will solve the determinis- 
tic equation for the first moment whereas for higher order 
moments we will take into account the lowest order cor- 
rection due to noise. 

Let us consider a process with cubic drift [l3| also 
called Bernoulli oscillator, 

dx = {-ax - bx 3 )dt + V2CdW t , (32) 

where b, C > 0. When a < the system represents a 
stochastic motion in a bistable potential. The statisti- 
cal moments of this stochastic system were early studied 
within a variational approach context [l4j . However, this 
method lands up introducing an extra (fitting) parame- 
ter which we adamantly want to avoid, since it would 
increase the level of uncertainty of the results. 

Figure [3] exhibits the first, second and fourth finite- 
time KM moments, obtained for different values of the 
sampling interval r. Notice that, with increasing r, D\ 
tends to a linear form that can erroneously lead to assume 
that the drift is linear. However, under this assumption, 
the inversion of D\ should give incongruous results for 
Dx, thus leading to discard the linearity of D\. For small 
noise intensity, the deterministic expression fitted (with 
two fitting parameters a and b) to the observed D\ allows 
to recover D\ . Notice also the complex influence of Di on 
the higher order coefficients, leading to forms with two 
minima, in contrast with the linear drift case. 

By integration of Eq. ([8]). one obtains the deterministic 
expression, 

D^M = ( -l)^, (33) 

that does not preserve the original simple cubic form un- 
derscoring the unique character of the linear drift. In the 
limit b — > 0, one recovers the linear case studied above. 
The case where b is a perturbative parameter may be 
of interest for systems in the vicinity of a phase transi- 
tion, such as in Refs. [l5{. In the limit a — > 0, Eq. (|3"3")l 
becomes, 

D**(xo,t) = ( , L 2 - 1 p . (34) 

Independent of the signal of a, for short r, Df e '(xo,r) is 
equal to —ax — b x 3 . For large values of t, the picture 
depends on whether a is positive or negative. In the 
former case, which corresponds to a single well potential, 
Dl det is defined by a straight line, 

Df et (x ,T) = -^, (35) 

T 

which corresponds to the independence limit. Mean- 
while, in the latter case, which entails bi-stable poten- 
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FIG. 3: Color online. Impact of r on the first, second and 
fourth finite-time KM coefficients, for the process defined by 
Eq. ([32]), with a = -0.1, b = 1 and C = 0.0025. The values of 
t are indicated on the figure. Symbols correspond to numer- 
ical estimates. Black solid lines correspond to the intrinsic 
coefficients Di = — ax — bx 3 , D2 = C and D4 = 0, while black 
dashed lines to the independence limit forms at r = 200. 



tials, one obtains, 



Df(x ,r) 



3 |o|t 



v/1 + fo^g(e 2 l-l - — l)/|a| 



(36) 

In this case, for large values of \a\ r, one can describe the 
limits corresponding to large values of b/\a\ x^e 2 a T , 



£>H*o,t) 



(37) 
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FIG. 4: Color online. Impact of the noise intensity on the first 
and second finite-time KM coefficients, for the same process 
of Fig. In all cases r = 5, and the different values C used 
are indicated on the figure. Black solid lines correspond to the 
intrinsic coefficients, black dashed lines to the independence 
limit at t = 5, and black dotted lines to the deterministic 
forms (J9J) at r = 5. Symbols correspond to numerical esti- 
mates and the associated colored lines to the approximations 
for small noise intensity given by Eqs. (|42[) and (|45[) . respec- 
tively. In the lower panel the inset shows the same data in 
log-linear scales. For the smallest noise intensity two different 
initial conditions were used to fill both potential wells. 



and for small values of the same quantity, 



£>r(zo,r) 



exp [|o| t] — 1 6exp[3|a|r] 



-x 



2 lair 



4- (38) 



As a result, we can verify that the polynomial depen- 
dence of Df et (xo , t) is preserved for the central region 
that dwindles as the sampling rate increases resulting in 
the straight line limit (|3T|) . which describes the full de- 
pendence (beyond relaxation scale) situation. 

In the generic case D\ = —hx n (with h,n > 0), the 
solution to Eq. (|5]) is, 



D**{x ,t) 



1 



(1 + (n 




(39) 



which includes Eq. (|34p as a particular case (for n = 3) 
and also the linear case with 6 = (for n = 1). In 
the latter instance it is provable that the exponential 
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functional dependence is recovered. For all the cases and 
in the limit r — > 0, one recovers Di(x) whereas in the 
opposite limit r — > oo. one gets —x/r. 

For small fluctuations, higher order finite-T KM mo- 
ments are dominated by the deterministic part given by 
Eq. ([9]) with further corrections dependent on the noise. 

For the cubic drift, in the presence of small fluctua- 
tions, such that C ~ a, 



(x) = (x) det +ci + 0(C 2 ,aC,a 2 ) 



with, 



ci = —CbxoT 



2 3 + 4y + V 
(1 + lyfl 2 



(40) 



(41) 



and y = bx^r. Then, 



D 1 =Df et + — + 0{C 2 1 aC,a 2 ), (42) 

T 

where D\ et is given by Eq. (|33|) . 

Similarly, for the second moment we found the first 
correction to the deterministic part 

(x 2 ) = ((x) det ) 2 +c 2 + 0(C 2 ,aC, a 2 ) , (43) 



with 



Then, finally 



C2 



2Ct 



(l + 2yy 



(44) 



= l(Df et ) 2 + C2 ^ 0Cl +O(C 2 > aC,a 2 ), (45) 



2t 



which verifies that D 2 tends to D 2 = C as r — >• 0. 

Figure 0] shows the impact of the intensity of noise on 
the finite-time estimates. In contrast to the linear drift 
case, here the noise influences Dx, in such a way that 
as the intensity of noise C increases, D\ approaches the 
independence limit even for small sampling rates. This 
could be attributed to the decorrelating role of noise (it 
is worth regarding that D\ already depends on the noise 
for non-linear potentials). Similarly to the linear case, 
D 2 largely disagrees with the true form D 2 , departing 
towards the independence limit, for large noise intensity, 
and towards the deterministic limit for small noise. No- 
tice that if the intensity of noise is too small, one of the 
wells of the potential may become inaccessible. Precisely, 
let us focus on the insets in the lower panel of fig. [4] 
When one is in the no-fluctuation regime the system is 
characterized by three fixed points, namely one unstable 
(at x = 0, which prevails if and only if x = is the ini- 
tial condition) and two stable points (at x = ±y/a/b). 
Minimal noise is able to make the observable, which we 
can envisage as a virtual particle, move off the unsta- 
ble point to one of the wells, each one defined by one 
of the remaining extremes (minima), but it gives a very 
little rate of transition between the regions of stability. 



As the noise intensity augments, the rate of transition 
increases as well and we are able to detect the finite val- 
ues of the KM moments in the vicinity of the unstable 
point. This can be adduced by computing the first non- 
vanishing eigenvalue, A, when the FPE is transformed 
into a Schrocdinger equation ■ In first approximation, 
the eigenvalue A of this bistable potential, is given by, 



r°° 2 
A = / dxe^ x + ~^ x 



dye 



(46) 



For the values presented in Fig. [H A goes from A w 0.2448 
(for c = 0.1) to A w 6 x 10~ 13 (for c = 10~ 4 ). 



V. FINAL REMARKS 

In this manuscript we obtained analytical results about 
the impact of the sampling rate on the KM coefficients 
directly computed from timcserics. 

We analyzed stochastic processes subject to a linear 
drift, which already comprise a huge variety of processes. 
We managed to compute exact expressions for the evo- 
lution of conditional moments, from which we obtained 
the finite-time KM coefficients and correlation functions. 
The first moment and the linear correlation function are 
both independent of the noise intensity. 

We also analyzed standard nonlinear drift cases, where 
the scenery exhibits a sharp change. The moments 
are now ruled by a cascade of differential equations 
which hampers the obtention of exact analytical solu- 
tions. Nonetheless, approximate expressions valid for 
small noise intensity were achieved by means of the ad- 
joint operator approach. Differently from the linear case, 
the intensity of the noise affects every conditional mo- 
ment and the correlation function as well. We also 
showed that by increasing the noise intensity the func- 
tional dependence of the KM heads towards the indepen- 
dence regime, which is characterized by a polynomial of 
the same order. This fact has notorious consequences in 
time series analysis, since one can be brought onto a sit- 
uation where independent-like coefficients are obtained, 
but one is unable to assign this finding to either an im- 
proper sample rating or a strong noise intensity. 

For sampling rates smaller than the relaxation time, 
in a steady state approach, our results permit one to 
sieve the set of processes giving raise to the heuristic 
stationary probability density function in order to select 
the appropriate Markovian differential stochastic process 
and consequently to determine the respective parameter 
values. Furthermore, it enables the judgement of the 
validity of the Markovian proposal by supplying precise 
estimates of higher order KM moments. Otherwise, if 
the relaxation time is shorter than the sampling rate, 
the measured KM coefficients will tend to a polynomial of 
order matching the order of the moment. This functional 
form is independent of the primary stochastic process, 
hence introducing uncertainty in the recognition of the 
process. 
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Uncertainties might be hedged by inspecting probabil- 
ity densities and correlation functions as well as other 
measures of dependence such as the relative (Kullback- 
Leibler) entropy. 
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